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ABSTRACT 



Context. Polar corona is often explored to find the energy source for the acceleration of the fast solar wind. Earlier observations show 
omni-presence of quasi-periodic disturbances, traveling outward, which is believed to be caused by the ubiquitous presence of outward 
propagating waves. These waves, mostly of compressional type, might provide the additional momentum and heat required for the 
fast solar wind acceleration. It has been conjectured that these disturbances are not due to waves but high speed plasma outflows, 
which are difficult to distinguish using the current available techniques. 

>A/Vt?s. With the unprecedented high spatial and temporal resolution of AIA/SDO, we search for these quasi-periodic disturbances 

in both plume and interplume regions of the polar corona. We investigate their nature of propagation and search for a plausible 

interpretation. We also aim to study their multi-thermal nature by using three different coronal passbands of AIA. 

Methods. We chose several clean plume and interplume structures and studied the time evolution of specific channels by making 

artificial slits along them. Taking the average across the slits, space-time maps are constructed and then filtration techniques are 

applied to amplify the low-amplitude oscillations. To suppress the effect of fainter jets, we chose wider slits than usual. 

Results. In almost all the locations chosen, in both plume and interplume regions we find the presence of propagating quasi-periodic 

disturbances, of periodicities ranging from 10-30 min. These are clearly seen in two channels and in a few cases out to very large 

distances 250'') off-limb, almost to the edge of the AIA field of view. The propagation speeds are in the range of 100-170 km s~^ 

The average speeds are different for different passbands and higher in interplume regions. 

Conclusions. Propagating disturbances are observed, even after removing the effects of jets and are insensitive to changes in slit 
width. This indicates that a coherent mechanism is involved. In addition, the observed propagation speed varies between the different 
passpands, implying that these quasi-periodic intensity disturbances are possibly due to magneto- acoustic waves. The propagation 
speeds in interplume region are higher than in the plume region. 

Key words. Sun: corona - Sun: oscillations - Sun: UV radiation - Sun: transition region - waves 



1. Introduction 

During solar minimum, polar corona is structured with sev- 
eral bundles of open magnetic field regions called plumes. The 
dark regions between plumes are called interplumes. Plumes 
and interplumes are characterized by slightly diflTerent proper- 
ties. Plumes are denser and cooler than interplume regions (e.g., 
[ IWilhelml l2006l) . I nterplume regions exh ib it broader line pro- 
files than plumes (iBaneriee et al.l l2000bh . lOfman et all (Il997h 
first reported the quasi-periodic variations in polar plumes us- 
ing the white light channel of UVCS/SoHQ They proposed that 
they ar e signatures of compressional waves. Deforest & GurmanI 
(Il998l) identified the ubiquitous prese nce of these propagating 
disturbances using EIT/SoHO. OfmanetaLl &99, 2000) in- 
terpreted these disturbances as slow magneto-acoustic waves. 
Similar work has been performed by m any others using spec- 
troscopi c observations of CDS/SoHO (iBanerjee et al.l l20QQal 
12001 alibi: lO'Shea et al.^ "2006", "2007') and more recently using 
SUMER and EIS ( Banerjee et al. 2009; Gupta et al. 2010), and 
all concluded that these disturban ces were magneto-aco ustic 
waves. From STEREO observations jMcIntosh et al.l(l201Ql) sug- 
gested that these are quasi-periodic high speed outflows rather 
than waves. A clear understanding of the basic physical mecha- 



nism involved is essential to quantify the significance and con- 
tribution of these quasi-periodic disturbances to the acceleration 
of the fast solar wind and/or coronal heating It has also bee n de 
bated whether the plumes (^Casalbuoni et al.|[r999 : Gabriel et al 
2003, 20 01) or interplum es dwilhelm et al 
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l2QQ3i: iGupta et al.ll201Ql) are the preferred channels for the ac- 
celeration of the fast solar wind. In this letter, we use the data 
from AIA/SDO in three coronal channels, to explore and study 
propagating quasi-periodic disturbances using space-time maps, 
constructed from several plume and interplume regions and try 
to address these problems. The simultaneous full-disk imaging 
capability of AIA/SDO in different channels with high spatial 
and temporal resolution, made this study possible. In the follow- 
ing sections, we first discuss the observational data used in this 
study, then the techniques used for the analysis and results ob- 
tained, and finally draw our conclusions. 



2. Observations 

We use the imaging data fro m the Atmospheric Im aging 
Assembly on-board SDO (AIA: iBoemer. P.. et al.ll201Ql) in the 
Extreme Ultraviolet regime in three coronal channels centered 
at wavelengths 171 A, 193 A, and 211 A. These three coro- 
nal channels correspond to their emissions mainly from Fe ix, 
Fexii, and Fexiv and their temperature responses peak at 0.8 
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MK, 1.25 MK, and 1.6 MK, respectively. Hereafter, we use the 
central wavelengths to refer to these channels. Data used in cur- 
rent study is spread over two hours on July 20, 2010 from 02:00 
UT to 03:59 UT. It is Rice compressed, calibrated level 1.0 data 
and is directly used in this analysis without applying any further 
corrections. We did not perform any de-rotation since the rota- 
tion effect would be insignificant on our analysis. The AIA pixel 
size is ^0.6'' and the cadence is 12 s. Fewer than ten frames out 
of 600 frames from these two hours, in each passband are miss- 
ing, which are added by linear interpolation. We chose to study 
the south polar region because it contains clear plume structures. 



3. Analysis & results 

In the south polar region, several structures are identified 
following peaks (plumes) and depressions (interplume) in 
intensity using contours (see Fig.[T]). At these locations, artificial 
slits of 60 pixels (^^36'') width in plume region and 30 pixels 
(^W) width in interplume region are extracted from each time 
frame, averaged over width and stacked to construct space-time 
maps. We chose a total of 12 slits, seven in plume and five in 
interplume regions, covering the on-disk part of the plumes 
wherever visible and extended almost to the end of the AIA 
field of view. The intensity contours and chosen slits with their 
exact widths are overplotted with solid (plume) and dotted 
(interplume) lines on a snapshot in 171 A in Fig. [T] The slit 
numbers are marked at the bottom right of each slit. Diff'erent 
widths are chosen to isolate interplume regions from the plume 
regions, hence study them separately. Same locations are used 
in all the three passbands as shifts of a few pixels between 
them (because of any misalignment), if at all present, can be 
ignored. The space-time maps thus obtained for each slit are 
then processed by smoothing over time using 200 points (^40 
min) smoothing, subtracting from the original, and normalizing 
to the smoothed image. The smooth subtraction enhances the 
faint variations and the normalization will take care of the 
radial intensity fall. Any quasi-periodic disturbance in intensity, 
propagating along the slit will appear as alternate bright and 
dark ridges in this processed map. These processed maps for slit 
2 (plume) and slit 9 (interplume), can be seen in Fig.O One can 
see from Fig. [2] first panel, that one such disturbance traveled 
almost 250''off'-limb, reaching the end of slit, and might have 
traveled yet further. The slope of these slanted ridges indicates 
the projected propagation speed of the disturbance, and the 
spacing between two successive ridges implies the periodicity. 
Since these ridges are fainter and rarely seen in 21 1 A, we limit 
our analysis to 171 A and 193 A . The speeds are calculated 
by fitting straight lines to the ridges and periodicity calculations 
are performed using a wavelet analysis. The horizontal lines 
marked on Fig. |2] enclose the region averaged and used in our 
wavelet analy sis. We use the Morlet fu nction as the mother 
wavelet (see iTorrence & Compol Il998l) . The results of the 
wavelet analysis for slits 2 and 9 can be seen in Fig. [3] The 
top panels of each subfigure in Fig. [3] shows the original light 
curve which sometimes contains longer periods that may be 
due to background variations. These are removed by choosing 
the appropriate background from the running average (over 
50-60 min) of the original, which are overplotted. The Middle 
panel shows this background trend subtracted light curve. The 
description of other panels can be found in corresponding 
figures. The primary and secondary peaks in periodicities, 
repres ented as PI and P2 are w ell above the 99 % confidence 
level (iTorrence & Compoll 19981) for a white noise process. The 



Table 1. Observed periods and projected propagation speeds for 
all the slits marked in Fig.[I] in both 171 A and 193 A. 





Fe IX 


171 A 


Fe XII 


193 A 


Slit No. 


P1(P2) 


Speed 


P1(P2) 


Speed 




(min) 


(km s-^ 


(min) 


(km s-^ 


1 


24.3 (13.2) 


103.2+1.4 


28.8 (12.1) 


133.1+2.3 


2 


18.7(12.1) 


113.8+1.4 


20.4(17.1) 


140.2+2.3 


3 


31.4(13.2) 


110.9+1.5 


26.4 (13.2) 


124.2+1.6^ 


4 


18.7(11.1) 


117.5+1.9 


22.2 (10.2) 


142.2+3.2^ 


5 


28.8 (14.4) 


109.9+1.6 


26.4 (10.2) 


137.6+2.5 


6 


26.4(12.1) 


116.5+1.7 


24.3 (17.1) 


132.6+1.9 


7 


28.8 (18.7) 


121.5+2.1 


31.4(17.1) 


136.9+3.4 


jd±cr 




113.4+5.6 




135.2+5.6 


±6 




+0.6 




+0.9 


8 


24.3 (14.4) 


137.4+2.3 


24.3 (14.4) 


149.8+2.6^ 


9 


26.4(17.1) 


142.5+1.8^ 


26.4 (18.7) 


164.7+2.9^ 


10 


18.7(12.1) 


144.9+2.7^ 


18.7(11.1) 


153.9+3.6^ 


11 


24.3 (...) 


127.2+1.6^ 


22.2 (...) 


167.4+3.8^ 


12 


20.4 (13.2) 


121.2+2.1 


28.8 (15.7) 


165.1+3.9 


li±cr 




134.6+9.1 




160.2+7.0 


±6 




+0.9 




+ 1.5 



Notes. PI and P2 are the primary and secondary periods (see Fig. |3). 
11 is the average of all the values above, cr is the standard deviation 
representing statistical error, and 6 is the instrumental uncertainty prop- 
agated from individual values. Jet-like features that can be seen from 
the movies. Enhanced intensity regions showing a decelerating ridge 
at the end of space-time map. This is due to some eruption from a region 
above slits 4 and 10. 



periods and propagation speeds for all the slits are tabulated in 
Table [T] If the second peaks are not available or below the 99 % 
significance level, they are marked as The periodicities 

are in the range of 10-30 min and the propagation speeds are 
100-170 km s"^ The errors in the calculation of speed, are 
due to instrumental uncertainties, given by = v(y + 7), 
where v, x, and t are the propagation speed, and the vertical 
and horizontal extents of the lines marked in space-time maps, 
and in addition, 6x, and 6t are the pixel scale (0.6'0 and time 
cadence (0.2 min), respectively. The errors are small, owing to 
the high spatial and temporal resolutions of the instrument, but 
the true uncertainty in the slope measurement might be larger 
than this. We may alternatively consider the width of the strip, 
in time direction as 6t, but despite choosing the narrower ones, 
the uncertainties are too large and the actual human error would 
be definitely smaller than this. Hence, we consider the standard 
deviation (cr) over the average propagation speed (yu) of all the 
values as the measurement error. Considering this, there is a 
noticeable diff'erence in the average speeds between plume and 
interplume regions and also between diff'erent passbands. 

Mostly in the 193 A passband, we detected in a few cases, 
one or two ridges that are steeper and brighter than other ridges. 
In one instance, we also observed a decelerating bright ridge at 
the end of the time sequence, which is brightest in 211 A. The 
former case is identified as the eff'ect of strong jets and the latter 
is due to some eruption that occurred just above the slits 4 and 
10 and propagated off'-limb. The slits aff'ected by these events 
are marked in Table [T] The eruption can be seen most clearly 
from the processed movies available onlineQ. Strong jets should 
be identifiable in the original movies. We emphasize original 
here, since the movies, processed to extract the fainter variations 



^ [ftp : //ftp ■ iiap . res . in/krishna/aia_polar] 



2 



S. Krishna Prasad et aL: Propagating disturbances in polar corona 



/yA 171A 2D-July-2010 02:00 UT 




Fig. 1. South polar region of the Sun as seen 
by AIA/SDO through one of its EUV chan- 
nels centered at 171 A. This snapshot is taken 
at 02:00 UT on 20 July 2010. Overplotted dot- 
ted curves are the contours of three different in- 
tensity levels indicating clear plume and inter- 
plume structures. The rectangular boxes delin- 
eate the locations of artificial slits extracted for 
the analysis, following intensity contours. Slits 
over plume regions are 60 pixels (^36'') wide 
and those over interplume regions are 30 pix- 
els (^18'') wide. These are marked diff'erently 
using solid lines for plume locations and dot- 
ted lines for interplume locations. Slit numbers 
are labeled at the bottom right of each slit. The 
movies show the temporal evolution as seen in 
the 171 A, 193 A, and 211 A channels. 



Processed space-time map in 171 -slit 2 Processed space-time map in 193-slit 2 Processed space-time map in 171 -slit 9 Processed space-time map in 193-slit 9 




Fig. 2. Space-time maps with time on X-axis, constructed from slits 2 and 9 (See Fig.[T]) and processed, in two coronal channels of 
AIA, 171 A and 193 A. First two panels are for slit 2 (plume) and the last two are for slit 9 (interplume). The slanted solid line, 
marked in each panel, following the ridges is used for the propagation speed estimation. The horizontal dotted lines in each of these, 
enclose the rows averaged for wavelet analysis. Corresponding maps in 21 1 A are shown in Fig. [4] (available online only). 



(< 5%) may produce the same visual impression for both waves 
and jet outflows. We note that the fainter jets, if at all present, do 
not affect our analysis because of our wider slits. 



4. Discussion 

We discuss various possible explanations of our results and a 
plausible interpretation. We chose wider slits than usual to avoid 
the efl'ects of fainter (few % above background) jets, if any were 
present at all. It is appropriate here to note that, though we 
present the case of 60 (plume) and 30 (interplume) pixel widths, 
we found the same results for slits as wide as 90 pixels (^^54'') 
and as narrow as 15 pixels (^9'0 which implies that a coher- 
ent mechanism is involved. Since the jets cannot be produced so 
coherently, the wider slits will average out the efl'ect of fainter 
jets, but the strong jets may remain. Using simple calculations, 
we roughly estimate that a jet of a five pixel (~3'0 width and 
of intensity 10 % above background will become a < 1% vari- 
ation assuming there is one jet at a time inside the 60 pixel slit. 
Hence, the stronger (> 10%) jets may appear in our space-time 
maps, but can be easily identified from the movies and filtered 
out. We found few such examples, most of which were eas- 
ily identified from the presence of distinct ridges with steeper 
slopes that are indicative of higher speeds. We also came across 
a faster jet that was indistinguishable from the other ridges be- 
cause of its inclination to our slit 9. These cases increase the am- 
biguity and shall be avoided by a careful inspection of movies. 
Following these steps carefully, we avoid the case of jets, yet still 
find propagating quasi-periodic disturbances in almost all cases 
that more importantly are insensitive to changes in slit width 



within certain limits. [Deforest et al.l (1 19971) suggested that the 
ionization temperature of the plume and interplume regions is 
around 1.0-1.5 MK. The propagating disturbances we find can 
also be seen clearly in 171 A and 193 A but rarely in 211 A, 
whereas the jets are seen mostly in 193 A and 211 A . This 
might be because jets are streams of hotter (than background) 
material whereas the (compressive) waves are local modulations 
in density. Another point is, though the periodicity values range 
from 10 min to 30 min there being three periodicities in most 
cases, with the first one varying between ^12 min and ^^18 min, 
the second one in-between 18 min, 24 min, and 30 min, and 
the third one around 45 min (though, inside the cone of influ- 
ence, see global wavelet plots in Fig.O. The repeated presence 
of these selective periodicities might be indicative of harmon- 
ics. Furthermore, the diflference in the propagation speeds be- 
tween 171 A and 193 A cannot be explained by a multi-thermal 
jet or outflow scenario whereas the temperature dependence of 
acoustic speed can explain the diflTerence. The acoustic speeds 
are related to temperature by the formula, Cs ~ 152 T^^^ m 
with T in K, as given by iPriestI (Il984l) . The ratio of observed 
speeds is 1.19, which can be compared to the acoustic speed ra- 
tio 1.25 considering the peak temperatures of 171 A (0.8 MK) 
and 193 A (1.25 MK). Now, comparing theoretically estimated 
acoustic speeds in 171 A (136 km s"^) and 193 A (170 km 
s"^) to those observed average values listed in Table [T] it is 
suggested that these disturbances are of slow magneto-acoustic 
type. However, if the actual uncertainties in speed calculations 
are larger (about 20%) then the diflTerence in speed will be neg- 
ligible and the interpretation in terms of acoustic type may not 
be valid. We also attempted to use the 211 A and 171 A pair. 
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Fe IX 171 Light curve — Original 



Fe IX 171 Light curve - Original 




Fe IX 171 Light curve — After trend subtraction 
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Fe IX 171 Light curve - After trend subtraction 
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Fig. 3. Wavelet analysis results for slits 2 and 9 (see Fig.[T|) in 171 A. Top two panels in each of these show the light curves, original 
and background-trend- subtracted, for the region enclosed by horizontal lines in Fig.O The bottom left panel shows the wavelet plot 
with contours enclosing the 99 % confidence regions, for a white-noise process. Bottom right panel is the global wavelet plot with 
99 % global confidence level, overplotted as a dotted line. The periods of primary and secondary peaks are also written in the text 
(middle-right). Corresponding plots in 193 A are shown in Fig. [5] (available online only) 



since their acoustic speed ratio is large. But from the space time 
maps corresponding to 21 1 A line, the ridges are clearly visible 
only for slit 2 location (see Fig |4]). The measured propagation 
speed for slit 2 in 21 1 A is 142.3+4.4 km s"^ and the speed ra- 
tio with respect to 171 A is 1.25 compared to the theoretical 
value 1.41. For the plume and interplume comparison, the prop- 
agation speeds are slightly higher in interplumes but we found 
no noticeable diff'erence in the ridge pattern. In some interplume 
cases, though there are slight curvatures at the bottom, close to 
the limb, they could be more due to projection effects t han ac- 
celeration unlike the case reported by lOupta et al.l (l2010h . 



5. Conclusion 

We have traced several plume and interplume structures us- 
ing contrast-enhanced images recorded in three coronal chan- 
nels 171 A, 193 A, and 211 A of AIA/SDO . The excellent 
spatial and temporal resolution of AIA/SDO, using space-time 
maps, has enabled us to find ubiquitous presence of quasi- 
periodic disturbances, which are unaltered by any changes in 
the slit width within certain limits. This implies that these are 
coherent phenomena. We have suppressed fainter jets using 
wide slits and excluded the effects of stronger jets by analyz- 
ing the movies. We emphasize that although our observations 
imply that these disturbances can be mostly due to presence of 
waves, the flow scenario cannot be completely ruled out based 
on these imaging observations. To explore the exact nature of 
these propagating disturbances, one needs to use coordinated 
imaging and spectroscopic ob serv ations, as demons trated by 
iDe Pontieu & Mclntoshl (l2QlQh and iTian et af] j201lh . We find 
that these disturbances travel far off-limb (^^250'0 and may pro- 
vide additional momentum to aid the acceleration of fast solar 
wind in coronal holes. We do not find any evidence of either ac- 
celeration or deceleration in either plume or interplume regions, 
but the propagation speeds are higher and close to the acoustic 
speed in interplume regions. 
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Processed space-time map in 211 -slit 2 Processed space-time map in 211— slit 9 

300 

250 




200 



^ 100 




150 



Time in nnln. starting at 02:00:14.05 on 20-July-10 



Time in mln. starting at 02:00:14.05 on 20-July-10 



Fig. 4. Processed space-time maps in 21 1 A, for slits 2 (left) and 9 (right). The ridges are faintly seen in slit 2 up to ^90'' off-limb. 
The slanted solid line overplotted on the left map, following the ridge, is used for propagation speed calculation. 
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Fig. 5. Wavelet plots for slits 2 (left) and 9 (right) in 193 A. The description of individual panels is given in Fig. [3] 



